Chemical Potential and the Nature of the Dark Energy: The case of phantom 
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The influence of a possible non zero chemical potential fi on the nature of dark energy is inves- 
tigated by assuming that the dark energy is a relativistic perfect simple fluid obeying the equation 
of state (EoS), p = up {ui < 0, constant). The entropy condition, S > 0, implies that the possible 
values of to are heavily dependent on the magnitude, as well as on the sign of the chemical potential. 
For i-L > 0, the cij-parameter must be greater than -1 (vacuum is forbidden) while for fj. < not only 
the vacuum but even a phantomlike behavior {lj < —1) is allowed. In any case, the ratio between 
the chemical potential and temperature remains constant, that is, fj,/T = fio/To. Assuming that the 
dark energy constituents have either a bosonic or fermionic nature, the general form of the spectrum 
is also proposed. For bosons /i is always negative and the extended Wien's law allows only a dark 
component with lj < —1/2 which includes vacuum and the phantomlike cases. The same happens 
in the fermionic branch for /x < 0. However, fermionic particles with ^ > are permmited only if 
— 1<LU< —1/2. The thermodynamics and statistical arguments constrain the EoS parameter to 
be a; < —1/2, a result surprisingly close to the maximal value required to accelerate a FRW type 
universe dominated by matter and dark energy (w < —10/21). 

PACS numbers; 98.80.-k, 95.36.+X 
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I. INTRODUCTION 

The current idea of an accelerating Universe is based 
on a large convergence of independent observations, and 
its explanation constitutes one of the greatest challenges 
for our current understanding of fundamental physics 
[l|, . In the context of a Friedmann- Robertson- Walker 
(FRW) cosmology, dominated by pressureless matter 
with density pm plus an extra component of density p 
and pressure p, the scale factor evolution is governed by 
the equation 3d/ a = —inG^pm + P + 3p). This means 
that a hypothetical component with large negative pres- 
sure may drive the evolution of an expanding accelerating 
Universe. This exotic component is usually termed dark 
energy or quintessence, and it represents about 70% of 
the total content in the Universe. The origin and the 
nature of dark energy is still a mystery, however, there is 
no doubt that its existence is beyond the domain of the 
standard model of particle physics Q. 

Among a number of possibilities to describe the dark 
energy component, the simplest and most theoretically 
appealing way is by means of a cosmological constant A, 
which acts on the FRW equations as an isotropic and 
homogeneous source with a constant equation of state 
parameter p/p = —1. On the other hand, although cos- 
mological scenarios with a A term might explain most of 
the current astronomical observations, from the theoret- 
ical viewpoint they are plagued with some fundamental 
problems thereby stimulating the search for alternative 
dark energy models driven by different candidates 0, Q ■ 

In the XCDM scenario, the dark energy component 
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is phenomenologically described by an equation of state 
p = up. The case lo = —1 reduces to the cosmologi- 
cal constant. The imposition a; > — 1 is physically mo- 
tivated by the classical fluid description [6j. This hy- 
phothesis introduces a strong bias in the w-parameter 
determination from observational data. In order to cir- 
cunvent such a difficulty, superquintessence or phantom 
dark energy cosmologies have been recently considered 
where such a condition is relaxed Q- In contrast to the 
usual quintessence model, a decoupled phantom compo- 
nent presents an anomalous evolutionary behavior. For 
instance, the existence of future curvature singularities, a 
growth of the energy density with the expansion, or even 
the possibility of a rip-off of the structure of matter at 
all scales are theoretically expected. Although possess- 
ing such strange features, the phantom behavior is the- 
oretically allowed by some kinetically scalar field driven 
cosmology [8], as well as, by brane world models and, 
perhaps, more important to the present work, a Phan- 
tomCDM cosmology is not ruled out b y th e present type 
la Supernovae and other observations [10, \lM ■ 

In a series of papers [H, have studied some 

thermodynamics and statistical properties of dark energy 
with no chemical potential {p, — 0). By using standard 
thermodynamics for a relativistic simple ffuid, we con- 
cluded that the case of phantom energy is ruled out be- 
cause the total entropy of a dark component with uj < —1 
is negative. In addition, by combining thermodynamics 
and statistical arguments the EoS was restricted to the 
interval —1 < lo < —1/2 and a fermionic nature to the 
dark energy particles was favored. 

Later on, thermodynamics arguments in favor of the 
phantom hypothesis were put forward by Gonzalez-Dfaz 
and Sigiienza [lj| . They claimed that the temperature of 
a phantomlike fluid is always negative in order to keep its 
entropy positive definite (as statistically required). This 
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new viewpoint was justified by arguing that the scalar 
field representation of a phantom field has a negative ki- 
netic term (fp' which quantifies the translational kinetic 
energy of the associated fluid system, and, as such, its 
temperature (a measure of the average kinetic energy) 
should be negative. Although temptative to some de- 
gree, both approaches have been considered in the lit- 
erature (see, for instance, [H, [l^ and Refs. therein). 
More important to the present work, they share a com- 
mon property, namely, the chemical potential of the dark 
energy fluid was set to be zero from the very beginning. 

In this article we reanalyze the thermodynamic and 
statistical properties of the dark energy scenario by con- 
sidering the existence of a non-zero chemical potential. 
In order to clarify some subtleties present in the earlier 
results, we rederive the physical quantities in the pres- 
ence of ^ by adopting the full thermostatistic approach 
proposed in Refs. |12| ITsj. This means that all thermo- 
dynamic and statistical properties follow directly from 
the EoS plus the hypothesis that /i ^ 0. In particular, 
the temperatures of the dark energy fluids must be al- 
ways positive definite. This is an interesting aspect of the 
present work since there are many scalar field representa- 
tions for dark energy fiuids, however, the thermodynamic 
laws are independent to what happens at a microscopic 
level as long as the equation of state has been defined. 
As we shall see, one of the main consequences of a neg- 
ative chemical potential is that the phantom scenario is 
recovered without the need to appeal to negative tem- 
peratures. In addition, a bosonic or fermionic nature of 
the dark energy component now becomes possible. 

The paper is planned as follows. In section 2 we discuss 
the thermodynamic constraints when a chemical poten- 
tial is introduced. In section 3 we consider a statistical 
analysis by assuming that the dark energy particles are 
massless and can have either a bosonic or a fermionic 
nature. In the conclusion section, a joint analysis involv- 
ing the thermodynamic and statistical constraints on the 
EoS w-parameter is presented 



II. COSMOLOGY, DARK ENERGY AND 
THERMODYNAMICS 



and negative (dark energy). The cases to = 1/3, 1, and 
— 1 characterizes, respectively, the blackbody radiation, 
a stiff- fiuid and the vacuum state while oj < —\ stands 
to a phantomlike behavior. 

Following standard lines, the equilibrium thermody- 
namic states of a relativistic simple fiuid are character- 
ized by an energy momentum tensor T"^, a particle cur- 
rent N°' and an entropy current S*" which satisfy the 
following relations 



T 



■Q/3 _ 



(p -I- p)u°'u 



a B 



■P9 



a/3 



S ',a — 0; 



(3) 



(4) 



(5) 



where (;) means covariant derivative, n is the particle 
number density, and s is the entropy density. In the 
FRW background, the above conservation laws can be 
rewritten as (a dot means comoving time derivative) 



p + 3(1 + Lo)p- = 0, fi + 3n- =0, s + 3s- ^ 0, (6) 
a a a 

whose solutions can be written as: 



|^ao-j3(i+w) 



P = Po\ — } , n = no , s = So , (7) 



where the positive constants po, no, so a-nd gq are the 
present day values of the corresponding quantities (here- 
after present day quantities will be labeled by the index 
"0" ) . On the other hand, the quantities p, p, n and s are 
related to the temperature T by the Gibbs law 



nTd{ — ) ~ dp — ^^-—^dn, 

n n 



(8) 



and from Gibbs-Duhem relation 18] there are only two 
independent thermodynamic variables, say, n and T. 
Now, by assuming that p — p{T, n) and p = p{T, n) and 
combining the thermodynamic identity [17| 



Let us now consider that the Universe is described 
by the homogeneous and isotropic Friedmann-Robertson- 
Walker (FRW) geometry (c = 1) 



ds^ = df--a^{t) 



dr"^ 



1 — nr' 



+ r^de^ +r^sin^ed(j}^ , (1) 



where k = 0, ±1 is the curvature parameter and a{t) 
is the scale factor. In what follows we consider that the 
matter content (or at least one of its noninteracting com- 
ponents is a fiuid described by the EoS 



p^ujp, 



(2) 



where p is the pressure, p is the energy density and uj a 
constant parameter which may be positive (white energy) 



(9) 



with the conservation laws as given by (6), one may show 
that the temperature satisfies 



— = f ^\ - = -3cj- 
T \dpj,ji a 

Therefore, assuming that w / we obtain 



(10) 



n = no 



^ T = To 



flo 



-3uj 



(11) 



3 



The temperatures appearing in the above expressions are 
positive regardless of the value of io. In particular, in the 
radiation case (lo = 1/3), one finds aT = QqTq as should 
be expected. As compared to this case, the unique differ- 
ence is that the dark energy fluid (even in the phantom 
regime) becomes hotter in the course of the cosmological 
adiabatic expansion since its equation-of-state parameter 
is a negative quantity. A physical explanation for this be- 
havior is that thermodynamic work is being done on the 
system 

It should be stressed that the derivation of the tem- 
perature evolution law presented here is fully indepen- 
dent of the entropy function, as well as, of the chemical 
potential /i. The above expressions also imply that for 
any comoving volume of the fluid, the product T'^V re- 
mains constant in the course of expansion and must also 
characterize the equilibrium states (adiabatic expansion) 
regardless of the value of /i. Further, by inserting the 
temperature law into the energy conservation law ([7]), 
one obtains the energy density as function of the tem- 
perature 
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FIG. 1: The allowed intervals of a; values (heavy lines) and 
forbidden (dashed lines) for null, positive and negative chem- 
ical potentials. Note that a large portion of the dark branch 
< is always thermodynamically permitted. However, for 
fj, > 0, the phantomlike behavior (a; < — 1) is thermodynami- 
cally forbidden. 



Pa 



(12) 



Now, in order to determine the chemical potential and 
its influence on the thermodynamics of dark energy, we 
consider the Euler relation [l^ 



Ts = p 



(13) 



where /i in general can also be a function of T and n |19l . 
[20j . By combining the above expression with equations 
(HD, d?]) and (dl]) we obtain: 



P = P-a\ — 



T 



where 



Pa = — [(1 + ^)pa - TqSo] 
no 



(14) 



(15) 



This straightforward thermodynamic result has some 
interesting consequences. In principle, the chemical po- 
tential may be either positive or negative, and it also 
depends on the values of the EoS cj-parameter. In partic- 
ular, fi is always negative (^o < 0) in the case of phantom 
energy, and becomes even more negative in the course of 
time (T grows with the scale factor during the cosmic 
evolution). It is also known that /x is zero in the case 
of photons (uj = 1/3) because they are their own an- 
tiparticles [2l|]. In this case, (jlSp yields correctly that 
SsqTq = 4po as should be expected. In general, if /i = 0, 
necessarily the relation sqTq = {l + uj)po must be obeyed, 
which is just the present day expression of sT = (1 +U!)p 
as required by (fTT 



At this point, the fundamental question is: How the 
chemical potential modifies the entropy constraints [T^ . 
[T3I I derived in the previous papers? 

In order to show that we compute explicitly the en- 
tropy of dark energy for a comoving volume V. As 
remarked before, the entropy function should scale as 
S ocT^V. Actually, 



5(r, v) = sv 



(1 +uj)pa - pono 



l/u 



V = soVo, 
(16) 

which remains constant as expected (see discussion below 
Eq. lfTT]) ). However, in order to keep the entropy S > 
(as statistically required), the following constraint must 
be satisfied: 



UJ > tJ„ 



Po 



(17) 



which introduces a minimal value to the w-parameter, 
below which the entropy becomes negative. This is a re- 
markable expression and its consequences are apparent. 
For instance, consider that /^o = (no chemical poten- 
tial). In this case, the smallest value of the uj— parameter 
is LOmin = ~1 and the previous analysis by Lima and Al- 
caniz |l3| is fully recovered, that is, the phantom domain 
(lu < —1) is thermodynamically forbidden. However, for 
a negative chemical potential, the phantomlike regime 
becomes thermodynamically allowed thereby recovering 
the hypothesis of phantom energy without appealingto 
negative temperature as proposed in the literature [14| . 
Note also that for a positive chemical potential not even 
a cosmological constant {uj = —1) is possible. In figure 
1, we summarize the basic thermodynamic results. 
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FIG. 2: Solutions for the bosonic case (jj, < 0). a) For given values of the pair (x,y), the straight lines and curves represent 
the l.h.s. side and the r.h.s. of (|2ip . respectively. The intersection points between the curves and the straight lines represent 
the desired physical solutions for ui > 0. The standard radiation solution (/j, = and cu = 1/3) is indicated by a black 
point, b) Curves representing the right hand side of (|2ip for different u values in the region u) < —1/2 (lower half plane) and 
— 1/2 < Lj < (upper half plane). The cosmological constant case {lj — —1) is indicated. For both diagrams the fugacity 
y = exp(^/fcflT) < 1. 



III. DARK ENERGY WITH A CHEMICAL 
POTENTIAL: STATISTICAL BEHAVIOR 

Another interesting feature of a dark energy compo- 
nent with a non zero chemical potential is related to its 
spectral distribution. The generalized Wien-type spec- 
trum for dark energy with U = has already been dis- 
cussed in the literature [12, [3 ■ A different approach for 
the phantomhke behavior involving the modulus of the 
temperature has also been proposed In the present 
case, since the temperature is positive, we simply add 
the chemical potential /i to the spectrum previously de- 
rived [H, • More precisely, we postulate the following 
spectral distribution: 



ye 



■2oj 



(20) 



where x = hc/ksXT and y — exp(fi/kBT) = 
exp(/io/^B21)) is a constant fugacity. When /ip = the 
above expressions reduces to the one obtained in p^ . 
The solution of the above algebraic equation can be de- 
rived both numerically and graphically. We are only in- 
terested in solutions with positive a;, because the tem- 
perature is always positive. Due to the presence of the 
chemical potential, the analysis of the above condition 
will be done separately for bosons and fermions. 



,\lu 



(18) 



where e = +1 stands for the Fcrmi-Dirac distribution 
and e = — 1 to the Bose-Einstein one, and a is a positive 
and w-dependent constant. Here it is important to note 
that for bosons the chemical potential is always negative 
or null, while for fermions it may be positive or negative 

0. 

A direct consequence of is related to the displace- 
ment Wien's law. The wavelength for which the distri- 
bution attains its maximum value is determined by the 
condition 



AniT' — 



he 



1.438 



where x'{uj) is the root of the transcendent equation 



(19) 



A. Bosons 

Let us now analyze the bosonic case (e = — 1). It 
proves convenient to rewrite condition (|20p in the follow- 
ing form: 



\ny — x = \n [ 1 — 



1 -I- 2uj 



(21) 



For each value of /i, the left hand side (l.h.s.) of the 
above expression is a collection of straight lines with slope 
equal to —1. Since /x is always negative or null for bosons, 
it follows that < y < 1 so that — oo < Inj/ < 0. This 
means that the l.h.s. of (PTjl is a collection of parallel 
straight lines on the lower half plane. The first straight 
line is the trivial solution with zero chemical potential. 

Note also that the right hand side (r.h.s.) of (PT|) in- 
volves the singularity for cu = —1/2, and, as such, must 
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FIG. 3: Curves representing the right hand side of (|22p . for 
different u values in the region uj < —1/2 and uj > 0. 



be separately analyzed. Apart from this point, we have 
3 different intervals, namely: u > 0, uj < —1/2 and 
-1/2 <LU <0. 

In Figures 2a and 2b we display the main results. 
When w > we have a collection of curves represented 
in Fig. 2a. All of them cross the straight lines in some 
point X > thereby indicating the solutions of the al- 
gebraic equation (|2ip . The standard radiation solution 
(/i = and = 1/3) is indicated, however, it should also 
be remarked the theoretical possibility of a radiation so- 
lution with /i 7^ 0. For uj < —1/2 we have the collection 
of curves represented on the lower half plane of Figure 2b, 
superposed to the straight lines. This analysis show that 
all these curves cross the straight lines in some positive 
X value, indicating a solution to the algebraic expression 
(I^TI) . Finally, on the interval — 1/2 < a; < we have the 
set of curves represented on the upper half plane of Fig- 
ure 2b thereby showing the absence of physical solutions. 

In summ, a simple graphic analysis shows that there 
are two intervals of uj for which the condition l|2ip has a 
solution, namely, uj > and uj < —1/2, while the inter- 
val — 1/2 < (jj < is statistically forbidden. Therefore, 
unlike the previous analyzes with = [H, [ll| , the EoS 
UJ < —1/2 for bosons now becomes possible when the 
chemical potential is negative. This includes the phan- 
tom dark energy as a physical possibility. 



B. Fermions 

The analysis of the fermionic case (e = +1) is similar 
to the bosonic one, but now the chemical potential can 
be either positive or negative. The condition (|2ip for the 
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FIG. 4: Thermodynamic and statistical constraints. The al- 
lowed (heavy lines) and forbidden (dashed lines) values of uj 
for the bosonic and fermionic cases with < and /i > 0. 
The phantom branch uj < —1 is excluded for the fermionic 
case with a; > 0. Note that the dark branch —1 < uj < —1/2 
for bosons is now possible. 



fermionic case reads 



In 2/ — a; = In — 1 



l + 2uj 



(22) 



The analysis on the l.h.s. of is similar to the 
bosonic case, the unique difference is that /io can be pos- 
itive. In this case, we have y > 1 and \ny > 0. As in the 
previous case, the discussion on the r.h.s. of (|^^ depends 
on the UJ values. For the cases ui < —1/2 and w > we 
have the curves represented in Figure 3. We see that all 
curves crossing the straight lines for some positive value 
of X yield a valid solution for ([22| . Note also that on the 
interval —1/2 < uj < 0, all the curves are in the nega- 
tive a;-axis (negative temperatures), and, therefore, none 
of them cross the straight lines (two reasons for the the 
interval be a forbidden region). 



IV. CONCLUDING REMARKS 

In this paper we have investigated the thermodynamic 
and statistical properties of a dark energy fluid with 
equation of state, p = ujp, by assuming that its chem- 
ical potential is different from zero. 

In Figure 4, we summarize the main results of the 
present analysis by combining both approaches. As dis- 
cussed in section 2, the regions with 5 < are thermo- 
dynamically forbidden. Note also that many dark energy 
fluids satisfy the combined constraints regardless of the 
fi sign, that is, a large interval of negative uj values is al- 
lowed by thermodynamic and statistical considerations. 
However, a phantom like behavior (w < — 1) is permitted 
only for /i < 0, and the corresponding massless particles 
can have either a bosonic or fermionic nature. 
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It was also proved (see also Fig. 4) that the EoS param- 
eter of a dark energy fluid obeying a generalized Wien's 
law always satisfies the constraint uj < —0.5 (a thermo- 
statistics limit). This upper limit is surprisingly close to 
the maximal value of the EoS w-parameter necessary to 
accelerate the present universe. Actually, in order to ac- 
celerate a FRW universe dominated by matter and dark 
energy, the EoS parameter must satisfies the inequality, 
UJ < —(1 -f- ^Im/^x)/^- Therefore, for il„i = 0.3 and 
il^ = 0.7, as indicated by the present observations 
one finds the dynamic constraint uj < —10/21. 

Finally, it should be stressed that for /i = one finds 
i^min = — 1 (see Eq. (|T7)) ) in accordance to the results 
previously derived by Lima and Alcaniz [13] . The present 
analysis with /i 7^ also opens the possibility for an EoS 
parameter w < — 1, thereby recovering the idea of a phan- 
tom dark energy without to appeal to negative temper- 



atures. Perhaps more interesting, unlike the results for 
/i = which favored only a fermionic nature to the dark 
energy fiuid (phantom and nonphantom), it was demon- 
strated that a bosonic kind of dark energy becomes pos- 
sible from a thermostatistics viewpoint. 
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